MKPH-T-00-08 

Analysis of the low-energy r/NN-dynamics within a three-body formalism * 



A. Fix a and H. Arenhovel fc 
a Tomsk Polytechnic University, 634034 Tomsk, Russia 
b Institut fur Kernphysik, Johannes Gutenberg- Universitdt Mainz, D- 5 5099 Mainz, Germany 

(February 8, 2008) 

The interaction of an 77-meson with two nucleons is studied within a three-body approach. The 
major features of the r/NN-system in the low-energy region are accounted for by using a .s-wave 
separable ansatz for the two-body r)N- and iViV-amplitudes. The calculation is confined to the 
(J 71 ; T) = (0~; 1) and (1~; 0) configurations which are assumed to be the most promising candidates 
for virtual or resonant TyiVJV-states. The eigenvalue three-body equation is continued analytically 
into the nonphysical sheets by contour deformation. The position of the poles of the three-body 
scattering matrix as a function of the ^^-interaction strength is investigated. The corresponding 
trajectory, starting on the physical sheet, moves around the r/NN three-body threshold and con- 
tinues away from the physical area giving rise to virtual ryiV TV-states. The search for poles on the 
nonphysical sheets adjacent directly to the upper rim of the real energy axis gives a negative result. 
Thus no low-lying s-wave ^TVJV-resonances were found. The possible influence of virtual poles on 
^ ' the low-energy jjiViV-scattering is discussed. 
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' I. INTRODUCTION 

The low-energy interaction of ?/-mesons with few-nucleon systems has recently become the subject of vigorous 
investigations. A great deal of attention has been devoted to the search for possible bound or resonance states in these 
systems. After an experimental study at Brookhaven yj], which casts some doubt upon the possibility of observing 
77-nuclei with relatively large mass numbers A>12, the attention was redirected to the interaction of 77-mesons with 
light nuclei (2[j^] . Recent theoretical investigations have mainly been stimulated by the precise measurements of the 
photoproduction and hadronic reactions with 77-mesons in final states. The obtained results are often interpreted as 
strong experimental hints that correlated states for light 77-nuclear systems might exist (see e.g. Q). The case in point 
is the strong energy dependence of the experimental cross sections observed near the 77-production threshold |p|,io|. 

In the present work we deal primarily with the dynamical properties of the r/N TV-system. The understanding of 
these properties is important for the following reasons. Firstly, the ^AW-system is interesting in itself. It is the 
simplest 77-nuclear system which admits an exact solution within the three-body formalism. Furthermore, it may be 
regarded as an example of a three-body system in which the pairwise driving forces are attractive and which may 
in principle develop a bound (or virtual) state, or a low-energy three-particle resonance near zero energy. Secondly, 
the study of r^AA-scattering may serve as a promising tool for the investigation of the corresponding processes in 
complex nuclei. The possible existence of such correlated objects would point to the crucial significance of two-nucleon 
mechanisms in the ^-interaction with nuclei. This fact in turn would require the revision of the simple first order 
approximation to the ?y-nuclear optical potential adopted by the currently available models 

The main features of the ^AA-dynamics in a quasideuteron state of spin, parity and isospin ( J ff ; T) = (1~; 0) are 
already elucidated in the literature. Ueda was the first to carry out an extensive calculation of ^-scattering within 
a three-body approach || and found that the rjd-systeni can form a quasi-bound state with a mass of 2430 MeV and 
a rather small width of about 10-20 MeV. Other calculations, also looking for bound or resonant ?yd-states within 
the Faddeev theory, have recently been reported by Shevchenko et al. Q as well as by Garcilazo and Pena ||. The 
authors of Ref. Q confirmed qualitatively the results of Q and fixed the values of the ^A-scattering length for which 
the previously bound ryd-state becomes a low-lying s-wave resonance. In Ref. [|| the possible existence of r]d bound 
states within the different r]N- and AA-interaction models is also studied. As for the experimental investigations, 
we would like to mention the measurements presented in |To| ] for the reaction np — > rjd where a visible increase of 
the ?7-meson yield over the mere phase-space calculation was observed near threshold. Recently, Metag et al. |9) have 
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reported new results for the cross section jd — > r/X. These authors also note a strong enhancement of events in 
the energy region of a few MeV above threshold. Such a feature cannot be explained within the truncated multiple 
scattering approach developed in Jl3], where only single r]N- and AW-rescatterings were included in the calculation 
as the most important corrections to the impulse approximation. This discrepancy is not surprising, since the validity 
of the latter approach is associated with the short range nature of the rjN- and AW-interactions in relation to the 
characteristic internucleon distance in the deuteron. In this situation, when thinking about the whole picture of the 
ijNN dynamics, one may assume that, when the attractive forces pull all particles together so that the two-body 
potentials overlap, qualitatively new features of the resulting 7/AW-interaction may be expected. 

In this paper we will study the question whether the dynamical properties of the ^AW-system allow the existence 
of a bound (or virtual) state or a three-body resonances in the low-energy region. Although this question has already 
been covered partially in the above mentioned works |^,^,^| we would like to reinvestigate it using a fundamentally 
different method, based on a search for poles of the scattering matrix in the complex energy plane. It possesses 
several important advantages over the conventional procedure of solving the on-shell Faddeev equation. Besides 
being more convenient for determining the exact position of S-matrix poles, the method provides additional insights 
into the source of appearance of these poles close to the physical domain. Furthermore, this approach allows us to 
investigate the (J g ;T) = (CP; 1) state which is not touched upon in the literature to the best of our knowledge. As 
was established in M], this configuration plays the dominant role in the final state of the 7c? — ► rjX reaction near 
threshold. Therefore, if one attributes the anomalous behaviour of the 7c? — > r/X cross section reported in |J to 
a near lying S-matrix pole, one has to search it in the (J n ;T) = (CP; 1) state. Broadly speaking, in view of the 
exclusive character of the measurements S , the marked enhancement of the near threshold yield can be assigned to 
the r/d final state. But the simple estimation shows that in order to explain the observed results one would need an 
enhancement factor of about 30 what in our opinion seems unlikely. In this context, large attention is focussed on the 
( J 71 "; T) = (CP; 1) configuration in this paper. 

Our principal tool is the three-body approach realized within the Alt-Grassberger-Sandhas formalism Jl5| ]. Since 
the search for the S'-matrix poles requires the analytical continuation of the dynamical equation into the nonphysical 
sheets of the Riemann surface, the analytical form for the driving two-body interactions is needed. For this reason, 
we use a simple separable potential of rank one and restrict our consideration to only s-waves in the two-body 
subsystems. This approach is well justified by the empirically established s-wave dominance of the low-energy r/N- 
and NN- scattering. There exist strong experimental and theoretical evidences that the low-energy ^A^-interaction 
is dominated by the formation of the £11 (1535) resonance. Analogously, the 1 S'o and 3 Si poles determine to a large 
extent the nucleon-nucleon low-energy interaction. Thus the s-wave isobars in the rjN and AW two-body channels 
are expected to be the main source of the ijNN-ioices. Therefore we hope that our separable model, though being 
quite simple, will reproduce the major features of the ?7AW-dynamics. 

In Sect. H we briefly describe the three-body scattering formalism pertinent to the present problem. Some details 



connected with the driving two-body interactions as well as the main calculational formulas are given in Sect. Ill 



In Sect. pV| the procedure of analytical continuation of the scattering equation into nonphysical energy sheets is 
described and the strategy for the search of the S'-matrix poles is presented. The discussion of our main results and 
the conclusions are presented in the last two sections. 



II. GENERAL FORMALISM 



In this section we will briefly review the properties of the three-body equation which will concern us in the remainder 
of this paper. Our starting point is the conventional three-body scattering theory in AGS form Jl5| . The three channels 
comprising two interacting particles and one spectator are labeled according to the number of the spectator. As already 
mentioned in the introduction, we take for each channel as the off-shell two-body t-matrix a rank-one ansatz 

t k (E) = \f k )r k (E)(f k \, (1) 

and we will call such an interacting pair in the following "isobar". Here \f k ) is the vertex function, and the isobar 
propagator T k is defined as 

r k (E) = 2* , (2) 

l- lk {h\Gf{E)\f k ) 

(2) 

where G k (E) denotes the free two-body Green's function with E as total cm. energy of the two-particle system. 
The separable form ([!]) enables one to reduce the three-body problem to a set of coupled effective two-body equations 
for the transition amplitudes X%j 
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Xij{W) = (l-5 l] )Z l] (W)+Y / (l-hk)Z tk (W)T k (E k )X kj (W), (i,j = 1,2,3), (3) 

fc=i 

where W denotes the total three-body energy while the isobar propagator r k depends explicitly on the invariant 
energy E k of the two-body subsystem, which is a function of W and the kinetic energy of the spectator. The energy 
dependent driving terms are 

Zu(W) = (fi\G{W)\fi) , (4) 

with G{W) being the free three-body Green's function. When evaluating equation (Q) in momentum space and making 
a partial wave decomposition, it becomes a set of one-dimensional integral equations of Fredholm type. It is known 
that for an inhomogeneous Fredholm equation to be singular, in which case the transition matrix Xij has a pole, it 
is necessary and sufficient for the corresponding homogeneous equation 

3 

F ij {W) = ^2(l-S ik )Z ik (W)T k {E k )F kj {W), {i,j = 1,2,3) (5) 
fc=i 

to have a nonzero solution for the same value of the parameter W. Thus the problem of searching for poles of the 
scattering matrix reduces to that of finding those values W for which the Fredholm determinant D(W) of (||) vanishes 

D(W) = 0. (6) 

We now turn to the 7?iV7V-system for which we adopt the labeling k = 1, 2 for a nucleon spectator and k = 3 for 
the r\ spectator, and furthermore change the channel notation as 

1 = 2 — > N* for the ryTV-isobar plus spectator nucleon , , . 

3 — > d for the iVW-isobar plus spectator meson. ^ ' 

The identity of the channels 1 and 2 reduces the 3x3 system (|^) to a 2x2 system which has the following operator 
form in our new notation 



Fdd = 2ZdN*TN*FN*d , (8) 

Fn-cL — Z^-dTdFdd + ZN*N*Tjy*FN»d , (9) 

FdN' = 2ZdN*Tjy*FN*N» , (10) 

Fn*w = Zff^dTdFdN* + Zn'N'Tn-Fn'N* ■ (H) 



A close examination of (|g) through ([n]) reveals that we have two decoupled sets of coupled equations sharing the 
identical integration kernel. Indeed, the coupled equations (||) and @ are transformed into (|l0| ) and (Q) by the 
substitutions Fdd — * F^n* and Fn*cI, Fn* m* . Therefore, it is sufficient to consider only one set for which we choose 
the second one. After inserting (|l^) into ([ll]) the former equation may be written in the closed form 

Fn*N* — {Zn*N* + 2Zjy*d1~dZdN') 1~N* Fjy*jy* . (12) 

It is diagrammatically presented in Fig. [j]. In analogy to a homogeneous Lippmann-Schwinger equation, we see that 
the driving term Zn* jv» plays the role of a meson exchange AA*-potcntial while the second term in brackets gives 
the mechanism associated with intermediate iVW-interaction where the 77 acts as a spectator. 



III. TWO-BODY INGREDIENTS 



Now we will specify the separable rjN and NN scattering matrices which determine the driving two-body forces 
in our model. Since we restrict the pair-interaction to s-waves only, the vertex functions have a simple structure 
determined by two parameters, the coupling constant g k and the cut-off [3 k 

(p\h) = fk (p) = 9k F k (p) with F k (p) = -g^-j , (13) 

Pi +p 2 

where p denotes the relative momentum of the interacting pair. 
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For the s-wave i-matrix of the NN — d isobar 

t d (p,p',E d ) = f d (p)r d (E d )f d (p') . 
the following parametrization has been used 

16™ 1 



9l 



a(3 d -2' ld 2M N 1 
where a is the iViV-scattering length. For the propagator r d one then obtains 



T d {E d ) = - 



1 



2M N 



16tt [i(3 d + y/M N (E d ~ 2M N ) / 
The N A^-interaction parameters were taken from the low-energy rop-scattering fit of Yamaguchi |16| 



f3 d = 1.4488 fin" 



5.378 fm for the 3 S'i-state, 
-23.690 fm for the 1 So-state. 



Analogously, we use for the r]N=N* i-matrix the following ansatz 

t N . (p, p\ E N * ) = f$ (p)r N * (E n , )f$ (p') . 
Here the iV*-propagator is given in the form 

t n , (E N * ) = (e n , - Mo -E v (E N *)-'Er,(E N .f) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



where Mo denotes the bare mass of the Sii(1535)-resonance. The resonance self energies associated with the couplings 
to the ttN and rjN channels are determined by 



£j (E N , ) 



1 

2^2 



/# 2 (P) 



p 2 dp 



E N . - E N (p) - Uj(p) + ie 2ujj(p) ' 



(20) 



where uij(p) = y m j +P 2 denotes the energy of meson tt j" and E^(p) = Mm + p 2 /2Mn is the nonrelativistic total 
nucleon energy. The meson- N* vertices are defined by (j = n, r/) 



f$(p)=9%lF$(p), with F#(p) = 



,0') rpU) 



GO, 



Pn» 



(21) 



The N* parameters were chosen in such a way that the main ratios of the hadronic decays of the Sii(1535) resonance 
are reproduced. We have taken 



g$ = 2.0, «,$ = 1.5, 0® = 6.51m- 1 , 0$ = 4.5 fin" 1 . 



The bare Six mass Mo was determined by the condition 

M + Me(T, n (M*) + E,(M*)) = M* , 



(22) 



(23) 



where M* = 1535 MeV is the mass of the dressed isobar. The choice (22) gives for the total and partial decay widths 
at the resonance position 



T = 150 MeV, 



r 



7) N 



jN 



0.5, 



(24) 



which is reasonably consistent with the values given in the 1998 Particle Data Group listings |l7j . 

For the actual evaluation of ( |l2| ) we use a natural representation in which the total angular momentum J and the 
total isospin T are diagonal. The r]NN wave function is expanded into a complete set of the following basis states 



4 



\q; JMj TM T ) = Kfcr,-) £ £ SMs {{(a^s^jSMs) \q; LM L ) Y LMl (q) . (25) 

LM L SM S 

Here oi and Tj denote the spin and isospin of the individual particles forming the isobar with the corresponding 
quantum numbers and £j. The partial waves are normalized as 

(q; LM L \q'; L'M' L ) = \s(q - q')5 L , L 5 M , ML . (26) 
q- 

In ( p5| ) we have already taken into account that the two-particle forces act only in I = states. The momentum <f 
characterizes the relative motion of the correlated (ij')-pair and a spectator particle. It is projected onto the partial 
waves specified by the total orbital momentum L with projection Mi,. Taking into account the parametrization 
adopted for the two-body vertices, the expansion (^5|) turns the operator equation ( |l2[ ) into a set of one-dimensional 
integral equations for the partial waves \q; LMl). In the following we consider only states with total orbital angular 
momentum L = 0. Then the two allowed configurations are (J W ;T) = (0~; 1) and (1~; 0). Consequently, only the 
TV TV-state contributes to the three-body state (0~; 1), and correspondingly the 3 5i TV TV-state to (1~; 0). 

In order to apply the representation (^H|) to Eq. (|l2|), we only need the spin-isospin recoupling coefficients between 
the states j and k for which one has 

{{((j i Oj)s k a k \SM s \{(o k o i )s, J o : j}S'M' s ) = S s > s§m' s m s \J 1 {2sj + l)(2s 4 + \)W {a l a k a j S; SjS k ) , (27) 
({(r,r J )( l T t }TM r |{(T t r,)i j T J }r'l4) = 5 t >t5 m > t m t yf&tj + 1)(2U + l)W{nT k T 3 T- t 3 t k ) , (28) 

where W(j\ J4 ; jbjs) denotes the standard Racah coefficient. 

Taking the actual quantum numbers of the participating particles, we obtain the following explicit form of the 
driving terms appearing in 1l2|) 



(p; JT\Z N * N *\p'; JT) = V v (p,p', W) + xKfe?' W) . (29) 

Here the spin-isospin coefficients are \ = —1/3 for (0 _ ; 1) and \ — 1 f° r (1 _ ; 0). The meson-exchange potential acting 
in the L = wave reads 

n m +/ F^l(p' + Tr M^p)F { '!(p+ jjM^p') 
lS7r { 2cj J (\p + p'\)[W-E N (p)-E N (p')-u }j (\p + p'\)j 



with (j = 7r, 77) . For simplicity we use the nonrelativistic relative meson- nuclcon momenta in the arguments of the 
rcgularization vertex form factors. The driving term associated with the intermediate TV TV-interaction has the form, 
using E d = W — u> v (q) — q 2 /AMm as the total cm. energy of the interacting TV TV-pair, 



(p; JT\Z N * d T d Z dN *\p'; JT) = V d (p,p',W) 

OO 

2 f q 2 dq 



it J 2bJ n {q) 


where the nucleon-exchange potential is given by 



V N (p, q, W) r d (w - - g 2 /4Mjv) V N (p', q, W) , (31) 



V N (p,P , W) = — / — ' \ ' d(p-p). (32) 

8?r J W-E N (jp) -w v (P ) -E N (\p + p ' ) 



The integrals in (30) and (32) can be evaluated analytically. The corresponding expressions, being rather cum- 
bersome, are listed in the appendix. Finally, we present the partial wave representation of our basic homogeneous 
equation (O) 

F N * N *(p,W) = ~ I V n , n ,( P ,p',W)t n ,(W - E N (p') - p' 2 /2M*) F N » N ,(p',W)p' 2 dp' , (33) 

7T . 
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where the notation 



V N * N .(p,p',W) = 2V d (p,p , ,W)+V v (p,p',W)+ X V 7r (p,p',W) 



(34) 



is introduced, and the invariant mass of the iV*-isobar in (B4J) was evaluated as E N , = W - E N {p') - p' 2 /2M*. 
Anticipating a later result, we note that the contribution from the pion exchange potential V v in ( p4| ) is practically 
insignificant, and almost the whole attraction in the rj N TV-system comes from the first two terms with approximately 
equal strengths. 

The method for searching the zeros of the Fredholm determinant D(W) is to approximate the integral in Eq. fl33| ) 
by a finite sum, transforming it into an ordinary matrix equation. Then (|6|) may be written as an algebraic equation 



det 



2 / 

% PjCjV N *N' {p t ,Pj,W)T N * W - E N (pj) 



2M* 



0. 



(35) 



Here Cj are the weights for the chosen quadrature (in the present calculation we have chosen the Gauss quadrature) . 



IV. STRUCTURE OF THE RIEMANN SURFACE AND CONTINUATION INTO NONPHYSICAL 

SHEETS 

Our main concern is to find the zeros of the Fredholm determinant in the complex energy plane resulting in poles 
of the scattering matrix. The structure of the manifold Riemann surface for the rjNN-system in (J W ;T) = (1 _ ;0) 
channel is presented in Fig. |[ Shown are the physical sheet II i and the nearest nonphysical sheets II2, II3 directly 
adjacent to IIi. The sheet II 1 has a conventional analytical structure. Namely, it may have poles corresponding to a 
possible formation of three-body bound states, and the unitarity right-hand cuts. The former stem from the Cauchy 
type integrals inherent in the integration kernel in (|33|). There are two main cuts which determine the structure of 
the Riemann surface for the three-body problem in the energy region of interest: 

(i) A two-particle cut, beginning at the two-body threshold W v d = 2Mn + m v + td where td < is the deuteron 
binding energy. The cut arises from the propagator Td(Ed) which has a pole when Ed — 2Mjy = £d- This pole is 
associated with the elastic ?7-deuteron scattering and gives the well known two-body contribution to the common three- 
body unitarity relations. The corresponding square-root branch point W = W v d develops the two-sheet structure 
{IIi, II2} typical for conventional two-body scattering. 

(ii) A three-body cut, starting at the three-body threshold W v nn = 2AlN+m v . This cut is induced by the two-body 
right-hand cut in Td corresponding to the d — > JVJV break-up and by the similar cut in r^r» associated with the decay 
N* —>■ r/N. The corresponding nonphysical sheet is denoted as II3 in Fig. |2[ Furthermore, the driving terms V„, V v , 
and Vjsr in (|3(]) through (^) contain the well known logarithmic singularities which are analogous to the dynamical 
(left-hand) singularities of the nucleon-nucleon OBE potential and correspond to the exchange of a real particle which 
becomes possible when W > W v nn ■ For some values of p and p' these singularities pinch the real axis producing the 
additional three-body cut beginning at the branch point W = W 71 mn- In the partial wave representation, this point 
is of logarithmic type and gives rise to an infinite number of sheets at W = W v nn- The structure of the Riemann 
surface associated with the logarithmic singularities is in itself of no physical interest and is therefore not presented 
in Fig. |. 

All sheets are connected as depicted in Fig. |[ The sheet ±I 2 may contain the poles which show-up as correlated 
two-body rjd states (virtual or resonant). The analogous states which may be observed in the three particle scattering 
processes are located on the sheet II3. In the absence of an actual three-body scattering experiment, these states may 
occur as final states in reactions with, e.g., deuteron break-up, such as 7c? — > r/np. Following the terminology accepted 
in the literature we call II2 and II3 as two-body and three-body sheets, respectively. 

The singularity structure outlined above arises from the free three-body propagators and the poles associated with 
the bound states of two-body subsystems. It is also presented in detail in Refs. |t8|-p2] in relation to the three- 
nucleon problem. The Riemann surface for the rjNN scattering matrix is formally more complicated due to the 
following reasons: 

(i) Since there is a rather strong coupling between the r]N and irN channels in the energy region of the 5n(1535) 
resonance, all sheets depicted in Fig. g have an additional cut beginning at irNN threshold^]. The typical structure 
of the Riemann surface associated with the r/N *-* ttN coupling is presented in Fig. a. 



1 Due to the spin-isospin selection rules the ird channel does not appear in the configurations (J n ; T) = (1 ; 0) and (0 ; 1) 
considered in this paper. Therefore no additional two-body cut starting at the -nd threshold is present. 
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(ii) The three-body sheet II3 has an additional square-root branch point associated with the quasi-two-body NN 



threshold. Treating the bare N* mass as noted in Sect. IV (see Eq. (|23|)), we keep the corresponding complex cut 



well away from the relevant energy region and, therefore, will ignore it in the following considerations. 

The structure of the Riemann sheet for the r/NN scattering matrix in the (J 71 "; T) = (CP; 1) state is different from 
that for (J n ; T) = (1~; 0). Namely, due to the virtual character of the 1 5'o pole, the corresponding two-body sheet 
II2 is "glued" not directly to III as previously but to the three-body sheet II3. 

From ordinary potential scattering theory it is known that the poles of the scattering matrix corresponding to 
virtual or resonant states are located in the nonphysical energy domain. The continuation into this area may be 
performed by going around the branch point at threshold in order to get the virtual pole or directly from the upper 
half complex energy plane by crossing the real positive axis to reach the resonance pole. This continuation may clearly 
be done if the analytical expression for the scattering matrix is available. Otherwise, one has to use the dynamical 
equation written in terms of the real energies in order to continue it into the nonphysical area. In the former case, 
several recipes may be used (for a review of the methods see, e.g. 

Turning to the three-body problem we would like to mention the method of analytical continuation based on 
a contour deformation technique. This method, developed in [ |19[|2C|1 , was previously applied to the three-neutron 
problem and later to the study of the ZPnn- and Ann- interaction ]23|]. In this paper we extend this technique to the 
niVTV-interaction including the inelastic 7riViV-channel. Firstly, we would like to review briefly the basic details of the 
method. Let us consider the following Fredholm equation 

00 

^' ^ = / 2mW-J-£-(r + pr F «> WW < (36) 



n 



where the function f(p', W) does not have any singularity in the relevant energy domain. For simplicity we consider 
three particles having equal masses m and use nonrelativistic kinematics. The energy W is the total kinetic energy. 
The physical region is determined in the usual fashion: dieW > 0, SraW = e — > +0. Because of the singularities of 
the kernel for real W one can not directly continue the equation down through the cut < W < 00. In this case, 
shifting the integration path for p' to the position C in Fig. || (the angle is arbitrarily taken to be 7r/4) we can cross 
the real energy axis and enter the area on the nonphysical sheet. Here the variable p is also taken on the contour 
C . The parabolic borderline separates the analytical domain from the forbidden area where the denominator in ( |36| ) 
may vanish for some values of p and p' . If one uses relativistic kinematics the permissible area slightly narrows. The 
continuation procedure is of course justified when the contour being deformed does not hit any singularities of the 
integration kernel. It must be noted that the method allows one to uncover only the certain part of the nonphysical 
sheet. However, varying suitably the contour parameters (a, 8), one can get the major part of the nonphysical domain. 
We note also that in the present calculation the procedure described above was applied directly to Eq. (|3^) , since all 
singularities of the kernel are contained in the Fredholm determinant D(W). 



A. Continuation into the lower half of the sheet IIi 



From the diagrams shown in Fig. |3| it is clear that the lower half-plane of the sheet Hi in Fig. g is at once the 
lower half of the three-body nonphysical sheet for the irNN channel reached from above by crossing the cut between 
irNN and r\d thresholds. This area is indicated as (II) in Fig. ||. Therefore, when continuing into this domain, 
the logarithmic singularities of the three-body propagator in V 7r (p,p', W) present an obstacle. A slight shift of the 
integration path (analogously to the pattern in Fig. 0) into the fourth quadrant allows one to reach the area shown 
in Fig. ||. Some calculational problems within this method may occur when the energy W approaches the real axis 
above the r/NN three-particle threshold since the analyticity domain degenerates into part of a straight line. In this 
limit one has to handle the numerics with greater accuracy. 



B. Continuation into the sheet II2 



In order to continue Eq. ( pq ) into the sheet H2, we adopt the following procedure. Entering into this sheet is 
accompanied by the movement of the pole po of the propagator into the fourth quadrant of p' . This requires the 
contour deformation as is shown on Fig. |^. In restoring the integration path to its original position, we pick up the 
new term in the potential Vd(p,p' , W) corresponding to the residue of the integrand in ( |3l| ) at q = po 
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V d (p,p', W) = IJ ^T^V N (p, q, W) T d (w - U> v {q) - q 2 /AM N ) V N (p',q, W) 
o n 

lu v (po) + 2M n 



where 

Po 



+ \e d \ - yjm* + AM N (W + \e d \ - Mn)}^ , (38) 



(compare with (|3l|)). In this manner we may cross the two-body cut in the interval W vd < W < W^nn and enter 
into the lower half-plane of the sheet II2 (transition 2 — » 2' in the notation of Fig. ||). The new domain of analyticity 
is shown in Fig. ^. 

Continuation into the upper half-plane of the sheet II2 (transition 1 — > 1') may be performed in a like manner. In 
this case, the pole po crosses the integration contour when moving from below into the first quadrant of the p'-plane. 



C. Continuation into the sheet II3 



For the continuation into the sheet II3 from the region above the segment W v nn < W < 00 of the real energy axis, 
we use the same technique as was outlined in the case A. In order to reach the domain below the rjNN threshold it 
is necessary to take 9 > 7r/4. In this case, however, one has to avoid the singularities on the imaginary p'-axis which 
come from the regularization form factors as well as from the factors of the type 2u>j (p) in the integrands in (|3C|). 
Therefore, the rotation angle 9 must be restricted to values less than n/2. Thereby, we also do not encounter any 
problem with the poles of r d and t^* . Choosing the mesh points from the deformed contour C, the three-body energy 
W can be taken to pass above the three-body threshold W„nn an d down through the cut into the lower half-plane 
of the sheet II3 as is shown in Fig. [?] (transition 4 — > 4' in Fig. |J). 

The analytical continuation into the upper half-plane of the sheet II3 is easily possible if the pion-exchange potential 
is ignored. Otherwise the logarithmic singularities from Vn bar the continuation path. For this reason, when entering 
this area from the lower half-plane of III, we will switch off the pion-exchange forces. In this case, in order to make the 
transition 3 — > 3' in the notation of Fig. ^, one just has to shift the contour into the first quadrant of the integration 
variable. Taking into account the smallness of the contribution from V v as denoted above, we believe that neglection 
of the pion-exchange potential does not spoil significantly the quality of our results. 



V. RESULTS AND DISCUSSION 



In order to obtain a better insight into the rjNN low-energy dynamics, we have investigated the trajectories drawn 
by the S'-matrix poles on the Riemann surface when varying the interaction parameters. In doing this, we are guided 
by the following consideration. We presume that the pole trajectory, on which one expects to find a resonance or 
virtual three-body state closest to threshold, will also contain the deepest bound state eigenvalue. In this regard, 
we artificially enhanced the strength of the rjN TV-forces until the first bound state (i.e., the pole on the sheet IT) 
appears. Then we have approached the actual physical situation by weakening the interaction strength to the value 



determined by the parameters given in Sect. III. Following this variation, the pole moves on the Riemann surface and 
finally arrives at several positions, developing in this way the bound (virtual) or resonant three-body state. Clearly, 
it may also appear on a Riemann sheet far removed from the physical region and thus will not influence the real 
scattering processes. 

Turning back to the rjNN -system, we take as a varying parameter the coupling constant g^jl . As its physical 
value we consider that presented in (p2j). We will start with the (1~;0) channel, the quasideuteron case. Before 
proceeding further, in order to sec how the pole may behave when approaching the threshold region, we will ignore 

(tt) (11) 

for the moment being the 7rA r A r -channcl by setting g N ,—0. Then we took gj^i—A and found the bound state pole at 

Wpd = W v nn — 8.08 MeV. As expected, the weakening of g^l implies the motion of the pole towards the threshold 

region. At g^jl—2.5 it overtakes the r\d two-body threshold and passes into the two-body sheet II2 producing a virtual 
r)d state. Further weakening of the 777V-attraction pulls the pole back along the negative real axis away from the 
physical region. We would like to note that this shape of the trajectory is what we can expect naively from ordinary 
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two-body potential scattering. Owing to the absence of the centrifugal barrier, the s-wave pole develops a virtual 
state and not a resonance when passing into the nonphysical sheet (here we do not touch upon exotic cases like the 
multitude of s-wave resonances in a deep square well p4|). We may expect to encounter the same pole behaviour in 
a three-body case since our physical situation is in effect the two-particle interaction with one being complex. 

Now including again the ttNN channel, we see that its effect is to shift the starting position of the pole downwards 
from the real axis and thus to transform the real rjNN bound state into the itNN resonance (or, what amounts to 
the same, the ijNN bound state with a finite lifetime). Then decreasing g^J; and allowing the eigenvalue W po i to pass 
the nNN threshold, we reach the situation illustrated in Fig. || where some pole trajectories corresponding to the 
different values of g^} are shown. One can see that the widths of the hypothetical nNN bound states are sufficiently 
less than that of the S l n(1535) resonance in this energy region. This is due partially to the fact that jyiViV-system 
spends a large fraction of time in the n(NN) state (not in N(N*)), which in this region has no direct coupling with 
any configuration in the continuum. This effect explains also the anomalous decrease of the nNN resonance width 
when approaching the ^-production threshold. As one can see from Fig. with increasing g^l the trajectories are 
shifted downwards and the point where the pole meets the real energy axis moves to the right. For the physical value 
(7^2 =1.5 this point is located at W ~ IF,, n at +1.51 MeV. The trajectory crosses the real axis above the three-body 
threshold and passes to the three-body sheet II 3 (dashed curve with open triangles in Fig. ^). If one go around the 
three-body threshold W^nn and enters the upper half plane of the two-body sheet II 2 (dotted line in Fig. ||), one 
finds the extension of the trajectory also into this sheet (solid curve with filled circles). We re call t hat we neglected 



the pion-exchange potential when calculating the trajectory passing into the sheet II3 (see Sect. JVC). For this reason 
the parts of the last two trajectories lying on the sheet 111 are not identical. We see however that the difference is 
not significant. The pole positions corresponding to g^l=2.0 are W po i = W v nn — (3.00 — z 13.67) MeV on the sheet 
IT3 and Wp i = W v d — (4.39 — i7.22) MeV on the sheet n 2 . 

After these findings we have carried out a careful search for poles on the lower half-plane of the sheets II2 and II3 . 
Of course, we were interested primarily in the energy area just beyond the real axis where possible resonant states 
may occur. Our conclusion is that no poles appear in the physically interesting domain at least for reasonable values 
of the interaction parameters. As for the presence of a possible resonance pole on the three-body sheet II3, one can 
assume that it may be found beyond the complex unitarity cut generated by the NN* break-up singularity in rjy* . 
In any case, this pole, if it exists, is far removed from the relevant energy region and will hardly influence the near 
threshold nNN interaction. Thus within our model, the ^-matrix poles nearest to the physical region are situated on 
the upper half plane of the sheets H2 and II3 and produce virtual rather than resonant states in the nNN -system.. 

Now we will consider the (J 77 ; T) — (0~; 1) channel. We make the same variation of the N* coupling constants and 
obtain the trajectories depicted in Fig. [| Switching off the irNN channel, we found the corresponding trajectory 
moving around the nNN threshold and ending in the third quadrant of the variable W — W v nn on the sheet II3. The 
similar pole behaviour was noted also for the £~wi and Ann s-wave configurations in Ref. p3|. Increasing g^j} and 

varying g^l from 4 to 2, we find the general picture to be qualitatively the same as in the quasideuteron case. One 
notes a visible shift of the trajectories towards the higher energies, which is most likely explained by the weakening of 
the W-interaction in the singlet 1 5'o state in relation to the more attractive 3 <Si state. The trajectory corresponding 
to g}^}=1.5 brings the final position of the pole to the point W po i — W v nn + (14.06 + H7.19) MeV. 

What are the physical consequences that we can extract from the results above? First, as we just have shown, 
it is unlikely that a low-energy resonance will be found in the (0~; 1) or (1~;0) channels. On the other hand, the 
virtual poles on the sheets II2 and II3 may influence the value of the scattering amplitudes through its proximity 
to the threshold region. Thus the next logical step in this direction would be to calculate, e.g., the 7/d-scattering 
near threshold. This problem, when treated within a three-body approach, is rather complex in itself and is beyond 
the scope of the present paper. The corresponding investigations for the elastic and inelastic ^-scattering will be 
presented elsewhere (see also |4|||)- Here we adopt the approximate formula for the s-wave elastic two-body scattering, 
when the amplitude has a virtual pole near zero energy (see e.g. p5[|) 

47T 

a " d = ' (39) 

where p denotes the cm. rjd momentum and p po i its pole position, determined as 



Pvoi = \ 2A ^ dW7? (Wpoi — Md — m v ) , ^smppoi < , (40) 
y M d + m n 

with Md being the deuteron mass. In Fig. |l0| we show the result obtained for W po i = W vd — (4.39 — i7.22) MeV. One 
sees a strong enhancement of the cross section as we approach the threshold energy, an effect which is typical when 
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the system possesses a weak virtual state. Therefore we conclude that in the region of a few MeV above threshold an 
anomalous resonance-like behaviour can in general be expected in the cross section of reactions with ryA^A^-channels 
in the final state. 

Finally we must note that our results are in contrast to those of Shevchenko et al. |D who claimed the existence 
of s-wave ^-resonances . They used the on-shell solution of the nonrelativistic three-body equation and found the 
counterclockwise rotation of the rid scattering amplitude in the Argand plot. The reason for this disagreement is not 
completely clear to us. On the other hand, we would like to mention the qualitative agreement of our results with 
those presented in || where a similar behaviour of the ^d-scattering cross section in the low-energy region was found. 



VI. SUMMARY AND CONCLUSION 



In order to answer the question whether there might exist a low-energy resonance or virtual r/NN-state, we have 
studied the typical pole trajectories W po i(g^jl) for the TV-system in L = states. Decreasing the ?7iV-interaction 

parameter g^i, the pole overtakes the three-body threshold and moves into the upper half plane of the sheets II2 
and II3 adjacent to the lower rim of the two- and three-body unitarity cuts on the physical sheet IT. Since these 
areas are not directly connected with the physical one, we have to conclude that no three-body resonances can occupy 
these trajectories. Moreover, we can assume that the very pattern of the pole trajectories forbids a low-energy rjNN 
resonance independent of the actual values of interaction parameters. We expect this feature to be valid in general for 
(L = 0)-r]NN configurations, at least within the model of a type presented here. On the other hand, our results point 
to a possible explanation of the strong enhancement of the 77-production cross section near threshold. This effect may 
be assigned not to a resonance but to virtual TyiV TV-states, indeed generated by the poles on the sheets ±I 2 and II3. 

We would like to emphasize that our conclusion is based essentially on the model of a pure s-wave interaction. 
The typical feature of an s-wave interaction known from the familiar two-body problem is to develop a virtual rather 
than a resonance state. In this context, it will be of interest to enrich the present study by including p-waves in the 
partial wave decomposition (^5|). It must, however, be kept in mind that p- waves are energetically far less favourable 
in a low-energy regime. Therefore, we think that even though a pole may come sufficiently close to threshold, quite a 
strong enhancement factor from a p-wave contribution would be needed in order to make a possible p-wave resonance 
visible in the cross section. 

Another aspect, left out in the present work, is the so-called direct TV TV-interaction which is not automatically 
accounted for in the conventional three-body approach. Due to the short-range character of this interaction, associated 
with its static nature, the corresponding terms are expected to give only a small contribution when compared to the 
full scattering amplitude (with respect to the role of these terms in ATV-interaction see e.g. J2(|). However, it should 
be remembered, that three-body low-energy dynamics is known to be rather sensitive to the model variation. Thus 
even a small perturbation may change some of the present results. 
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APPENDIX 



For the driving terms in (|30|) and fl32| ) one finds the following analytic expressions: 
(a) meson exchange potential (j — 7r, 77): 



pp' (af + c 2 )(«2 + c 2 ) 



{ 



In- 



2(a{-ai 



2 1nK 



1 — In r — 



)+(!<- 2)]} 



(Al) 



where 
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c = W-E N (p)-E N {p') 



l \2 



rrij + (p ± p ) 
Mjv + mj (j)2 



A/ 



TV 



J ' Mjv'' Mjv + to 



Mjy + fftj (j) 2 



M 



A'* 



7 2 

p , 



Mjv Mjv + raj 



and ^4^ = at ± i a± for k = 1, 2. If ai = a?2 = a, then with ^4^ = a ± i a± 



Vj(p,j/,W) 



9^0^ ( Mjv + m 



16tt V Mjv 



pp'(a 2 + c 2 ) 2 



I V c -a+ A7 J 



3a 2 



? c - 



In- 



A+ A~ 



4a 3 (A+ A~y 
(b) nucleon exchange potential: 



2 5Re 



16irpp' 



(tZ2 — ai)(a2 + c)(ai + c) 



. s , ai + PP' / , s , a 2 + pp' . c + pp' 

x < (a 2 + c) m - (ai + c) In - (ai — a 2 J hi 

ai — pp' a,2 — c — pp' 



where 



rP 2 , 



_ Mjv+m, / (7))2 , 2 \ m, 

2m, V^* +P ) + 2(M N + m v y 
a 2 = p 2 +p 2 +p' 2 /A, 
c = M N (W - E N {p) - w„(p') - Mjv) - (j? + p' 2 )/2 . 
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FIG. 1. Graphical representation of the homogeneous equation ( |l2| ) for the r/NN three-body eigenvalue problem. 




FIG. 2. Structure of the Riemann surface for the rjNN scattering matrix in the (1~ ; 0) channel in the near threshold region. 
Shown are the physical sheet IIi and the adjacent nonphysical sheets II2, and II3. Two vertical dashed lines pass through the 
two-body (W v d) and three-body (W^jvjv) thresholds. The numbers i and i! indicate the identified rims of the associated two- 
and three-body unitary cuts. 




FIG. 3. Relationship of the four sheets for a two-channel rjNN-irNN problem. The point moving around the rjNN threshold 
passes successively all four sheets denoted here by the roman numerals. 
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FIG. 4. Analytical continuation of the integral equation ( |36| ) within the contour deformation method. The left-hand panel 
shows the new integration path C, determined by the length of a finite straight line a and a rotation angle 9 (here 9 = 7r/4). 
On the right-hand panel, the new domain of analyticity in the lower half-plane of the nonphysical sheet is shown. The energy, 
corresponding to the corner point A, may be evaluated as Wa — — ia 2 jm. 
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FIG. 5. The energy domain of the sheet TIi available for the analytical continuation. The parameters of the integration 
path C having the form shown in Fig. are: a=200 MeV/c, 9 = 15°. 




FIG. 6. Continuation into the two-body sheet IT2. On the left-hand panel we show the pole po of the propagator Td moving 
across the integration path when the energy W passes through the two-body cut on the interval W v d < W < W v nn- The 
parameters of the new contour are the same as in Fig. ^[ The energy domain reached on the sheet II2 is shown on the right-hand 
panel. 
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FIG. 7. The same as in Figs. |E] and ^ for the sheet II3. The parameters of the shifted contour are a=200 MeV/c, 6 — 80° 
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FIG. 8. The pole trajectories of the r/NN scattering matrix for the (J" ; T) = (1~ 
to the different choices of the tvNN* coupling constant: gr^, =0 (sts 

The legend shows the values of the r/NN* coupling constant g^l being varied along the individual trajectory. The trajectory 
shown by the dashed curve is continued into the three-body sheet II3. The other trajectories pass into the two-body sheet II2 
(see the text). 
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FIG. 9. The pole trajectories of the r/NN scattering matrix for the (J"; T) = (0~; 1) state. Notation as in Fig. 
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FIG. 10. Elastic r;d cross section evaluated according to the approximate formula (pSS|). 
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